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PREFACE 


Recent  studies  of  mass-transfer  rates  associated  with  various  re¬ 
entry  environments  have  pointed  out  the  need  for  understanding  the 
flow-fie\d  phenomena  related  to  massive  injection  into  boundary  layers. 
Perhape  the  simplest  problem  of  the  type  that  can  be  examined,  without 
recourse  to  large-scale  computation,  concerns  the  effect  of  large  mass 
injection  on  the  constant-property  laminar  boundary  layer. 

In  this  Memorandum,  the  principal  goal  is  to  consider  the  effects 
of  pressure  gradient  and  large  mass-injection  rates  for  certain  bound- 
ary-layer  flows,  and  to  explore  the  nature  of  the  shear  layer  which  may 
join  the  inviscid  outer  potential  flow  to  the  inviscid  flow  dominated 
by  mass  transfer  near  a  wall.  In  addition,  the  similar  problem  in  nat¬ 
ural  convection  is  considered,  where  temperature  difference  plays  the 
role  of  the  driving  mechanism  for  flow  acceleration. 
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SUMMARY 


The  problem  of  the  constant-property  laminar  boundary- layer  flow 
with  large  mass-injection  rates  and  favorable  pressure  gradient  is  con¬ 
sidered.  Approximate  solutions  for  the  structure  of  blowing-induced 
shear  layers  are  obtained  for  flows  that  satisfy  the  requirements  of 
Falkner-Skan  similarity.  For  small  pressure  gradients  (0  <  1/2),  the 
asymptotic  structure  is  shown  to  consist  of  a  viscous  shear  layer  im¬ 
bedded  in  &n  inviscid  flow. 

A  preliminary  analysis  of  the  natural -convection  problem  with  large 
injection  is  performed,  and  the  way  in  which  the  ideas  used  in  the  Falk¬ 
ner-Skan  analysis  may  be  extended  to  this  problem  is  suggested. 
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LIST  OF  SYMBOLS 

injection  rate  at  surface,  -f(0) 
Falkner-Skan  stream  function 
f/C,  inner  variable 
1  -  s 


first  term  of  outer  expansion 
Prandtl  number 

* 

Reynolds  number,  — 

distance  between  surface  and  axis  of  symmetry  for  axisym 
metric  flow,  or  distance  between  surface  and  reference 
plane  which  is  parallel  to  body  force  in  two-dimensional 
flow 

temperature 

velocity  in  boundary  layer 
reference  velocity  (Eq.  (24)  or  (25)) 
velocity  at  edge  of  boundary  layer 
parabolic  cylinder  functions 


df/<m 


U/U 


normal  velocity  in  boundary  layer,  -y'Pr  Re 
inner  velocity  variable 
inner  distance  variable 

/  * 

normal  distance  in  boundary  layer,  vPr 

distance  in  direction  of  toc?y  force  (height) 
2 


u 

first  term  of  inner  solution 
Falkner-Skan  pressure-gradient  parameter 
volumetric  expansion  coefficient 


v 


y 


-x 


Falkner-Skan  similarity  variable 
stream  function  (Eq.  (36)) 
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I.  INTRODUCTION 


In  the  absence  of  pressure  gradient,  the  flat-plate  boundary  layer 
"blows"  off  the  surface  if  the  injection  of  mass  into  the  boundary  layer 
is  sufficiently  high.  This  "blow-off"  point  occurs  as  a  singularity  in 
the  solutions:  it  is  not  possible  to  find  solutions  of  the  boundary- 
layer  equations  for  blowing  rates  larger  than  some  critical  value. 

For  flows  with  favorable  pressure  gradient,  no  such  critical  "blow- 
off"  injection  rate  seems  to  exist.  For  a  certain  class  of  such  flows, 
the  Falkner-Skan  wedge  flows,  it  is  possible  to  obtain  solutions  to  the 
boundary- layer  equations  for  arbitrarily  large  blowing  rates. 

For  more  general  boundary-layer  flows,  with  large  mass  injection 
and  favorable  pressure  gradient,  the  solution  of  the  equations  of  mo¬ 
tion  in  the  neighborhood  of  the  surface  is  essentially  invlscld,  and 
a  simplified  Bernoulli  equation,  obtained  by  neglecting  the  component 
of  the  kinei/c  energy  associated  with  the  normal  velocity  component, 
is  an  approximate  solution  of  the  equations  of  motion. 

A  remarkably  simple  solution  to  the  problem  of  the  incompressible 

similar  laminar  boundary  layer  with  asymptotically***  large  injection 

(2) 

was  proposed  by  Pretsch  in  1944.  The  Falkner-Skan  equation  for  wedge 
flows  with  favorable  pressure  gradient  (3  >  0) 

fimn +  ffiFi +  f*1  ■  $  ■ 0 

was  considered  subject  to  the  boundary  conditions 

f-qCO)  -  0  f(0)  =  -C  f^OO  -  1  (2) 


A  more  complete  discussion  of  these  points,  as  well  as  a  discus¬ 
sion  of  the  Falkner-Skan  flows,  is  given  in  Ref.  1,  p.  243. 

** 

See  Ref.  1,  p.  347,  for  the  development  of  these  ideas. 

Asymptotic,  In  this  case,  refers  to  blowing  rates  which  are 
large  but  can  still  be  considered  within  the  context  of  boundary-layer 
theory . 
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where  C  New  variables  were  introduced, 

h(0  h  f (Tl)/C 

C  •»  VC 


(3) 


and  a  simplified  equation  was  obtained  by  neglecting  terms  of  order 
l/C  .  The  resulting  equation,  essentially  the  inviscid  version  of  the 
Falkner-Skan  equation  omits  the  viscous  term  f^  ,  and  consequently 
it  is  possible  to  satisfy  the  two  boundary  conditions  at  7)  "  0,  and 
the  one  at  7]  -* 00 ,  by  a  discontinuous  solution.  Pretsch's  solution  is 
valid  in  the  region  between  the  zero  streamline,  which  separates  in¬ 
jected  fluid  from  free-stream  fluid,  and  the  wall,  but  it  gives  no  in¬ 
formation  about  the  outer  flow  between  the  ?  to  streamline  and  the  po¬ 
tential  flow. 

Several  properties  of  Pretsch's  solution  are  worth  noting: 

1.  The  velocity  at  the  bounding  streamline  is  given  correctly, 

in  that  f_  =  U/U  -•  1  as  f  -♦  0.  For  nonsimilar  incompressible  boundary 
*1  ® 

layers,  with  large  injection,  this  is  not  necessarily  the  case.  As 
the  bounding  streamline  is  approached  from  the  wall  with  arbitrary  fav 
orable  pressure  gradient,  U 


inner 

U  (x  *  0)  ®0,  since  U  ~  x11  and  U 
e  e 

line. 


^/u^(x)  -  U^(x  »  0).  For  wedge  flows, 
inner  0ri  t*ie  bounding  stream- 


2.  For  values  of  @  >  1/2,  the  pounding  streamline  is  approached 
from  the  wall  with  zero  shear  stress:  for  3  =  1/2,  with  finite  shear; 
and  for  3  <  1/2,  with  infinite  shear. 

In  this  Memorandum  we  re-examine  the  problem  posed  by  Pretsch  and 
indicate  how  his  solution  may  be  matched  to  a  simple  viscous  solution 
which  is  valid  in  the  outer  flow,  and  which  is  necessary  to  smooth  out 
the  discontinuity  in  shear  resulting  when  3  <  1/2.  This  outer  solution 
is  shown  to  be  a  uniformly  valid  approximation  to  the  inner  flow  as 
well.  The  blowing  rate  C  enters  into  the  outer  solution  only  as  a 
multiplicative  constant.  The  outer  flow  consists  of  a  viscoi  3  shear 
layer  imbedded  in  a  large  inviscid  mass-transfer  layer  only  for  small 
values  of  3* 
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In  addition,  we  indicate  how  the  inviscid  approximation  of  Pretsch 
can  be  applied  to  free-convection  boundary  layers  with  large  mass-trans¬ 
fer  rates  to  obtain  simple  solutions  for  various  flow  quantities  in 
terms  of  a  simplified  Bernoulli  integral. 
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II.  ANALYSIS 

We  shall  find  it  convenient  to  work  with  modified  Von  Mises 
variables  u(f),  f,  which  are  appropriate  for  similar  flows: 


f_  s  u 

U 


f  a  f 


In  these  variables,  Eq.  (1)  becomes 


(4) 


u(uu£)£  +  fuu£  +  |3(l-u)  =  0 


(5) 


The  boundary  conditions  become 


u  =  0  f  =  -C 


u  -*  1  f  -* 


(6) 


2 

Letting  u  =  z,  we  obtain 

Vz(Zff)  +  fz£  +  2p(l  -  z)  =  0 
z(-C)  =  0 

Z  (co)  -*  1 

For  the  inner  region  near  the  wall,  a  new,  independent  variable,  f  s  f/C, 
is  introduced.  The  resulting  equation  and  boundary  conditions  are 

•~yz(z££)  +  fz~  +  2g(l  -  z)  =  0  (7b) 

c 

z(-l)  =  0 
z(»)  -*  1 


t 
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An  inner  solution  that  is  valid  for  fixed  ?  and  1/C  -*  0  is  con¬ 
sidered.  The  first  term  of  the  inner  solution,  zQ(f),  is  the  one  pro¬ 
posed  by  Pretsch  and  results  from  the  inviscid  equation  obtained  by 
neglecting  the  terms  of  0  {l/C2}  : 

fzo~  +  26(1  -  zo)  =  0  (8) 

The  no-slip  condition  is  retained: 

Z0(-D  -  0  (9) 

The  solution  of  Eq.  (8) ,  subject  to  the  boundary  condition  of  Eq.  (9) , 
is 


zo  =  1  -  <-f)2&  (10) 

This  solution  is  sketched  in  Fig.  1. 

For  values  of  f  in  the  neighborhood  of  0,  this  solution  has  the 
~  ~28-l 

property  that  f  ^  -*  f  ,  and  consequently  there  is  a  radical  differ¬ 
ence  in  the  shear  profile  as  8  is  less  than  or  greater  than  1/2,  as 
noted  earlier.  The  outer  solution  of  Eqs.  (8)  and  (9)  is  u  =  1  for 
all  values  of  f  >  0.  The  viscous  outer  solution,  which  smooths  out 
the  shear  profile,  joins  a  potential  flow  at  large  positive  values  of 
f  to  an  inner  rotational  flow  at  large  negative  values  of  f,  and  corre¬ 
sponds  to  a  ’Vorticity-interaction"  problem,  where  the  vorticity  arises 
from  the  total  pressure  variation  along  the  wall  which  is  convected 
along  the  streamlines  of  the  injected  fluid. 

The  solution  of  Eq.  (10)  cannot  really  be  continued  to  f  >  0,  At 
?  =  0,  z  reaches  the  value  at  infinity.  However,  the  resultant  pro- 
file  has  a  corner  at  zq  =  1,  f  =  0.  This  corner  is  removed  by  construct 
ing  an  expansion  that  is  valid  near  z  =  1,  which  is  matched  to  the  first 
term  of  the  inner  expansion.  The  matching  takes  place  as  f  -*  0-  in  zq, 
and  f  -»  -oo  in  the  outer  expansion.  Although  the  second  term  in  the 
inner  expansion  has  a  singularity  as  f  -*0-,  it  can  be  shown  that  in 
the  region  of  matching,  ics  contribution  is  negligible. 


0  0.5  1  1.5 


Transformed  normal  distance  variable,  /3rf/ C 


Fig.  1— Typical  inviscid  velocity  profiles 
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Thus  a  new  variable  g(f)  is  introduced: 

z(f)  »  1  -  g(f)  (11) 

where 

(a/1  '  g)gff  +  fgf  -  2pg  =  0  (12) 

g(oo)  -*  0 

g(-C)  =  1 

This  suggests  that  the  first  term  for  the  outer  expansion  of  g  is 

gouter'lcj  80<f)+-"  <13> 

The  resulting  equation  for  gQ(f)  i8  then 

®off  +  "  2^o  =  0  (14) 

gQ(oo)  -*  0 

gQ(f  -  -»)  -  (-f)2p 

The  solution  to  the  equation  and  boundary  conditions  of  Eq.  (14)  is 

g0  -  A  [exp  (^|-)]  [up(2B+  1/2,  f)  +  BVp(2|3  +  1/2,  £)]  (15) 

where  U^(2g  +1/2,  f)  and  V^(20  +  1/2,  f)  are  the  parabolic  cylinder 
functions  described  in  Ref.  3,*  and  A  and  B  are  constants  of  integra¬ 
tion. 


* 


In  the  notation  of  Ref.  3,  U 

P 


s  U  and  V 

P 


V. 
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We  note  from  Ref.  3  that 

ttV(2|3  +  1/2,  x)  »  T(1  +  2p)  sin  rr(2p  +  1/2)  *  Up(2p  +  1,  x)  +  Up(2p  +  1,  -x) 

From  the  asymptotic  expansions  for  U  and  V,  It  Is  clear  that  the  con¬ 
stant  B  must  be  zero,  since 


Vp(20  +  1/2,  x)  x2S  1  +  0  +  ••• 


as  x  -*  +  « 


and 


Up(2p  4-  1/2,  x)  ~  e  4  x-1"2B 


1  +  0 


(!0+  ••• 


as  x  -*  oo 


The  asymptotic  behavior  as  x  -»  -co 


e4  (-x)20  r" 


Up(20+1,  x)~  r(l  +  20) 


1  + 


-fr) 


+  •  •  • 


fixes  the  constant  A  to  be 


r.q  ±  ,2fl? 

2n 


The  one -term  outer  solution  Is  thus 


8o(f)  ~{cV  ^  2n2^  e  4  V20  +  1/2»  f)  (16) 


.ii 


The  shear  function  corresponding  to  this  outer  flow  Is 


2fTTn  s  -  r  r(1  X-28)  V2B  -  1/2,  f) 


£ 

4 


V2tT 


(17) 
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wherc  the  recurrence  relation 


(2$  +  1/2,  x)  -  (l/2)xUp(2g  +  1/2,  x)  +  Up(20  -  1/2,  x)  -  0 


has  been  used. 

Since  the  inner  solution,  Eq.  (10),  is  contained  within  the  outer 
solution,  Eq.  (16),  then 


g0(C, 


£) 


(0 


20 


m  +  2b> 


£ 

‘4 


Up(20  + 


1/2,  f) 


(18) 


is  a  uniformly  valid  one-term  expansion. 

Figure  2  indicates  the  shear -function  variation  for  various  values 
of  the  pres sure -gradient  parameter  6*  The  shear  function 


(C) 


26 


it 


is  equal  to 


(C)2B  2£^(n) 

in  che  usual  Falkner-Skan  variables. 

We  note  that  for  6  <  1/2  (in  this  case  6  =  0.1),  a  genuine  shear 

(2) 

layer  appears,  as  suggested  by  Pretsch's  inviscid  solution,  but  for 
larger  values  of  the  pres sure -gradient  parameter,  the  shear  is  a  maxi¬ 
mum  at  the  wall,  and  the  boundary  layer  remains  on  the  wall  rather  than 
being  "blown  off"  into  a  shear  layer. 


5-4-3-2-10  I  2  3  4  5 

Stream  function,  f 

Fig. 2— Shear  variation  in  wedge  flows  with  large  injection  rates 
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III.  FREE -CONVECTION  BOUNDARY  LAYERS  WITH  INJECTION 


The  inviscid  approximation  may  be  applied  to  free -convection  bound¬ 
ary  layers  if  injection  rates  are  sufficiently  large. 

(For  convenience,  we  will  use  the  notation  and  coordinate  system 
of  Ref,  4.) 

For  a  free-convection  boundary  layer  on  a  two-dimensional  or  axisym- 
metric  body,  the  equations  of  motion  in  dimensional  form  are 

uux  +  vuy  =  Wyy  +  60g(I  -  T„K/l  -  r'2  (19) 


uT  +  vT  =  T 
x  y  Pr  yy 


(rou)x  +  (rov)y  =  0 


(20) 

(21) 


where  j  ■  0  for  two-dimensional  flow  or  1  for  axisymmetric  flow.  The 
appropriate  reduced  variables  for  a  free-convection  boundary  layer 
are 


7  =  A A 

v  =  y/l 


Pr  Re  y 


Pr  Re  v 


(22) 

(23) 


where 


and  U  is  some  reference  velocity, 

U*  -  y^gATL  for  Pr  -  0(1)  (24) 


(25) 
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In  these  reduced  variables,  the  equations  of  motion  become 


uu  + 

X 


vu—  =  Pr 
y 


JU-—  + 

yy 


e0g(T 


V 


•A 


K>' 


(26) 


uTx  + 


(rou)x  +  <ro^y  -  0 


(27) 

(28) 


The  boundary  conditions  for  asymptotically  large  blowing 


are 


u(x,  0)  =  0 


u(x,  cd)  =  o 


v(x,  0)  =  Cv  (x) 
w 


(29) 


T(x,  0)  ■  T  (x) 
w 


T(x,  »)  =  T 


where  C  -*  eo,  and  v  (x)  0(1). 

W 

Our  discussion  centers  around  mass-injection  rates  which  are  much 
larger  than  l/y/vr  Re*  but  are  still  sufficiently  small  for  the  boundary- 
layer  equations  to  be  valid. 

New  variables  are  introduced: 


*  v 
v  =  — 

C 


*  z 
y  -c 


(30) 

(31) 


and  we  consider  solutions  to  Eqs.  (26)  to  (28),  subject  to  the  boundary 
conditions  of  Eqs.  (29)  valid  for  C  -  »,  which  are  of  the  form 

u  =  u(x,  y* ,  C)  =  uq(x,  y*)  +  ... 
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v  -  v  (x,  y  ,  C)  =  CvQ(x,  y  )  +  ... 


T  «  T*(x,  y* ,  C)  =  T*(x>  y*)  +  . . . 


•fc  •/(  iff  ^  ^ 

The  equations  for  uq(x,  y  ),  vq(x,  y  ),  Tq(x,  y  )  become 


*  5u  .  ou  /- 

uosf  +  vorl-^-i„)v^ 

&y 


Su 


<ro^2 


(32) 


-  *  * 

*  ®o  * 
u  r-2  +  v  — X  -  0 
o  Sx  o  .* 

Sy 


(33) 


Mrii  )  3(rV) 


o  o 


O  O' 


dx 


3y 


(34) 


with  boundary  conditions 


vQ(x,  0)  =  vw(x) 


uQ(x,  0)  -  0 


uQ(x,  oo)  =  0 


(35) 


T*(x,  0)  =  T*(x) 


T  (x,  »)  =  T 

o  o 


The  solution  of  these  equations  and  boundary  conditions  can  be  obtained 

* 

easily  by  the  introduction  of  Von  Mises  variables  x,  i|r  as  independent 
variables. 

Introducing 


_  * 
21- 
ax 


*  j 
-v  rj, 


* 

dii 

* 

ay 


*  j 

u  rJ 
o  o 


(36) 
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and  noting  that 


JL  _  ,  *  j  a  a  a  *  j  a 

8;'vV'  *  = 


results  in 


a_ 

ax 


=  P0g(T*(x,  **) 


(37) 


ax 


=  0 


(38) 


It  is  convenient  to  use  the  variable  x  ,  which  marks  the  position  on 

•jjj 

the  body  where  the  i|r  streamline  crosses  the  body  surface 


^(x*)vw(x*) 


(39) 


In  the  x,  x  coordinates, 

* 

u  *  x  *  / - s - 

Y  (x,  x  )  =  J  Pog(T(x  )  -  Tj/l  -  r^(Tl)  dTl  (40) 
x 


Tq(x,  x  ) 


T  (x  ) 
wx  ' 


(41) 


where  the  no-slip  condition  at  the  wall  for  u  has  been  applied.  The 

•||f  ^ 

solution  for  u  (x,  y  )  is  given  parametrically  in  terms  of  x,  x  ,  where 
*  ° 
y  is  determined  from 


r^(x)y*(x,  x*) 


ro<S>vw<?> 
uQ(x,  5) 


(42) 
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where 


j  */  \ 

VoW 


I 


.  *.  .  * 

x  r0(x  )vw(x  )  dx 


u0(x,  X  ) 


<M) 


Note  that  y  (x)  =  Cy  (x)  is  the  location  of  the  dividing  streamline. 
k  k°  °  k 

For  y  >  y  .  u  ■  0  and  T  =  T  .  The  dividing  streamline  a  surface 

O  O  O  co 

of  discontinuity  in  both  temperature  and  velocity  if  the  temperature 
at  x  =  0  is  not  T  .  If  the  temperature  at  x  <=  0  is  T  ,  then  there  may 
be  a  discontinuity  in  either  rhear  or  temperature  gradient  at  yQ(x). 

We  observe  that 


f  VT^d,  =  Z(x)  -  Z(x*)  (44) 

x 

k 

where  Z(x)  -  Z(x  )  is  the  difference  in  height  of  the  points  on  the 

k 

body  surface ,  x  and  x  . 

Equation  (40)  is  e  simplified  Bernoulli  equation  which  can  hr. 
written  as 


f-  (x,  x*)  =  Pg[T*(x*)  -  T^  [z(x)  -  Z (x* ) j  (45) 

For  a  vertical  flat  plate,  to  give  an  example,  Z(x)  =  x.  Then 


and 


2 

jn  (x>  x*)  “  eog[T*(x*)  -  tJl  •  (x  -  x*)  (46) 


Cx>  =  f 


v  (x*)  dx* 

*  * 
UQ(x,  X  ) 


(47) 
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If  T(x  )  =  T  ,  constant,  and 


w 


v  (x*)  =  V  , 
w'  7  w’ 


constant,  then 


_ \ _ v£ 

2  •  (B08)1/2(I„  -  I.)  2 


(48) 


Using  the  tame  ideas  which  we  have  introduced  in  the  calculation  of 
the  viscous  correction  for  forced  flow,  it  should  be  possible  to  ob¬ 
tain  the  viscous  correction  which  smooths  out  the  discontinuity  at 

*fc 

Yq(x)  for  certain  wall-temperature  distributions .  However,  it  appears 
that  the  problem  of  constant  wa’l  temperature  would  require  a  full  dis¬ 
cussion  of  a  thermally  driven  free  shear  layer  at  y  =  yQ. 

Since  it  is  not  at  all  clear  that  the  natural -convection  problem 
with  large  injection  rates  has  technical  applications  or  significance, 
we  have  included  the  inviscid  analysis  only  to  indicate  how  the  usual 
ideas  of  boundary-layer  theory  for  both  forced  and  natural  flow  can  be 
simplified  if  problems  involving  larger  mass-transfer  rates  are  con¬ 
sidered. 

The  principal  simplification  results  from  the  observation  that  the 
flow  in  the  neighborhood  of  the  boundary  is  inviscid,  so  that  the  bound¬ 
ary  layer,  which  occurs  at  the  wall  for  small  mass -transfer  rates,  can 
now  become  a  shear  layer.  For  the  forced-flow  case,  the  skin-friction 
parameter  is  determined  completely  by  the  inviscid  solution,  while  for 
the  natural- flow  case,  the  heat  transfer  to  the  boundary  is  zero. 
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